Transport reversal in a delayed feedback ratchet 
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Feedback flashing ratchets are thermal rectifiers that use information on the state of the system 
to operate the switching on and off of a periodic potential. They can induce directed transport even 
with symmetric potentials thanks to the asymmetry of the feedback protocol. We investigate here 
the dynamics of a feedback flashing ratchet when the asymmetry of the ratchet potential and of the 
feedback protocol favor transport in opposite directions. The introduction of a time delay in the 
control strategy allows one to nontrivially tune the relative relevance of the competing asymmetries 
leading to an interesting dynamics. We show that the competition between the asymmetries leads 
f~| ' to a current reversal for large delays. For small ensembles of particles current reversal appears as 

O ' the consequence of the emergence of an open-loop like dynamical regime, while for large ensembles 

uJ , of particles it can be understood as a consequence of the stabilization of quasiperiodic solutions. 

S ■ We also comment on the experimental feasibility of these feedback ratchets and their potential 

applications. 
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■ . I. INTRODUCTION 

Q ' Brownian motors or ratchets are spatially periodic systems that are able to induce direct transport rectifying thermal 

O , fluctuations. Two conditions are generally sufficient for the emergence of direct transport in these systems: breaking 

of thermal equilibrium and breaking of spatial inversion symmetry [l| . These systems permit one to get an insight 

fvj into non-equilibrium processes and are receiving increasing interest also due to their applications in nanotcchnology 

> and biology j],i, 13 . 

^^ [ Flashing ratchets are devices that rectify the motion of Brownian particles by subjecting them to a spatially periodic 

OO . potential that is alternatively switched on and off. Open-loop flashing ratchets operate without regard to the state of 
t/* I the system (open-loop control) implementing a periodic or random switching to rectify thermal fluctuations by taking 
advantage of the asymmetry of the potential [J, [a, Q . On the contrar\^ feedback ratchets (or closed- loop ratchets) use 
information on the particle distribution of the system to operate [3, [^, M, Ej \lM, H, \lM ' ^^^ ^^^ asymmetry of the 
feedback control protocol is able to induce a directed transport even for symmetric ratchet potentials. For instance, 
in the so-called maximization of the center-of-mass velocity protocol [3] the controller switches on the potential only 
if switching on would imply a positive displacement for the center-of-mass position (i.e., if the net force with the 
potential on would be positive). Feedback flashing ratchets have been recently suggested as a mechanism to explain 
the stepping motion of the two-headed kinesin |14L In another context, a feedback scheme has been used to perform 
control of chaotic trajectories in inertia ratchets [l5|. 
cd ' Feedback flas hing ratchets could be experimentally implemented monitoring the positions of a set of Brownian 

particles [la , [ITI . [lg | and subsequently using the information gathered to decide whether to switch on or off a ratchet 
potential according to a giving protocol. This experimental design will have to deal with a finite time lag between 
the collection of the information about the state of the system and the action because of the time interval needed for 
the measurement, transmission and processing of the informatio n fig. [20l . Time delays in the feedback also appear 
naturally in complex systems with self regulating mechanisms (see |2ll . |22| and references therein) . It is also remarkable 
for the ability of controlling chaos and improving coherence in excitable systems under delayed feedback [23l . |24| . The 
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feasibility of nanotechnological feedback flashin g ra tchet devices and their performance under the presence of a time 
delay has been analyzed very recently in Refs. [12, [lj|. In those works, and also in previous ones 0, S B [l3, hH . 
the two sources of spatial asymmetry involved, namely the feedback control and the shape of the potential, cooperate 
with the aim of maximizing the performance of the system. On the contrary, in this paper we investigate the effects 
of the competition between the potential asymmetry and the control asymmetry in a delayed feedback ratchet. 

We have observed a rich dynamics that includes transport reversal. The inversion of the current direction upon the 
variation of the system parameters is a well-known phenomenon in Brownian motors that can be produced by varying 
the characteristics of the non-equilibrium fluctuations 25, 26] or the parameters of the time-dependent perturbation 
that drives the system out of equilibrium [271 [28l . vM . Wk [31| . It also appears in other ratchet-like systems, such as 
deterministic inertial ratchets [32ll33|. The phenomenon of current reversal has great importance in particle separation 
devices [3J|, and in biology systems [3a]. In our present study current reversal is achieved just by varying the time 
delay of the system. 

We start below with the description of the collective flashing ratchet and the delayed feedback protocol that we 
consider. In the next section. Sec. Illli the evolution equations of the system are solved by Langevin dynamics 
simulations and the rich dynamics encountered (transport reversal, quasi-periodic modes of oscillation, multistability) 
is analyzed. We finally review and further discuss in Sec. llVl the implications of the results. 

II. MODEL 

The feedback ratchet that we consider consists of N Brownian particles at temperature T in a periodic potential 
V{x). The force acting on the particles is F{x) — —V'{x)^ where the prime denotes the spatial derivative. The state 
of this system is described by the positions Xi (t) of the particles satisfying the overdamped Langevin equations 

-fx^{t) = a{t)F{x,{t)) + C,{t); i=l,...,N, (1) 

where 7 is the friction coefficient (related to the diffusion coefhcient D through Einstein's relation D ~ kBT/j), ^i{t) 
are Gaussian white noises of zero mean and variance {^i{t)^j{t')) = 2^kBTSijS{t — t'), and a{t) stands for the action 
of the controller. The feedback policy uses the sign of the net force per particle, 
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as follows: The controller measures the sign of the net force and, after a time t, switches the potential on (a = 1) if 
the net force was positive or switches the potential off {a = 0) if the net force was negative. Therefore, the delayed 
control protocol considered is 

^(,)^fe(/(t-r)) ift>., ^3^ 

[ otherwise, 

with O the Heaviside function [0(x) = 1 if a; > 0, else 9(a;) — 0]. We have used a sawtooth potential of period L, i.e. 
V{x) = V{x + L), height Vb, and asymmetry parameter a: 

.r, N f— f ifO< f <a. 

The height Vq of the potential is the difference between the value of the potential at the minimum and at the maximum, 
while aL is the distance between the minimum and the maximum consecutive positions (Fig. [T]). Thus when a < 1/2 
both the asymmetry of the potential and the feedback protocol favor transport in the same direction, whereas when 
a > 1/2 there is a competition between them. We consider here the latter case. 

III. RESULTS 

We have performed numerical simulations of the Langevin equations ([T]) by using an Euler-Maruyama scheme [36| , 
which reveals different dynamics of the collective ratchet for different ensemble sizes when the delayed feedback is 
present. We distinguish between few particles (including iV = 1 as a particular case) and many particles. Following 
previous works [7|, we refer to the few particle case when the average long-time limit velocity of the center-of-mass. 
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FIG. 1: Sawtooth potential [Eq. (J4j] of height Vb = SfcsT and asymmetry parameter a — 2/3. Units; L = 1 and ksT — 1. 

(icm), is greater for the non-delayed feedback maximization protocol than for the optimal open- loop protocol, and 
many particle case otherwise. Typically the frontier between these regimes corresponds to A^ = 10^ — 10"^ particles 
for potential heights of the order of SfcsT or greater. 

Let us begin our analysis with the non-delayed (r = 0) feedback ratchet. 

A. Non-delayed feedback ratchet 



For one particle (N — 1) the exact expression for the average velocity can be obtained by solving a Fokker-Planck 
equation with the proper effective potential that includes the action of the controller. This expression has been derived 
in [7|, and is indeed valid for any asymmetry parameter < a < 1. It gives a positive flux that grows as DVQ/ikgTL) 
for small potential heights (Vb < ksT), and tends to the finite value 2D/{a?'L) for large potential heights (Vb ^ ksT). 
The fact that the flux goes to a constant value for large potential heights is a direct consequence of the overdamped 
nature of the ratchet. We remark that enlarging the value of the ratio Vb/(fc_Br) corresponds to effectively diminish 
the intensity of the white noise that accounts for thermal fluctuations. 

For the collective ratchet compounded of a few particles an approximation for the center-of-mass velocity can be 
obtained assuming a purely stochastic behavior (see Ref. [?[ for details). As the magnitude of the fluctuations of the 
force are of the order of the inverse of the square-root of the number of particles the stochastic approximation predicts 
the (icm) ~ i/vN decay observed in our simulations for any asymmetries. In fact this qualitative behavior remains 
valid for any number of particles (including many particles) provided the asymmetry is a > 1/2. In this latter case the 
potential asymmetry acts against the feedback protocol, which tries to favor positive currents, and then the potential 
is turned on in very small intervals of time as the controller rapidly switches it off. See Fig. [51 Thus the systems 
dynamics is effectively stochastic, contrary to the many particle case with cooperating asymmetries, where switches 
are slower and allow the system to have enough time to evolve in a quasideterministic way [7[. Therefore when the 
potential asymmetry competes against the feedback the flux decays with the number of particles as I/VjV even for 
many particles, contrary to the much slower 1/lnA^ dependence observed when both asymmetries cooperate l7l|. 

In any case the non-delayed protocol always gives a positive flux because it only switches on when it implies a 
positive displacement of the center-of-mass position. 



B. Delayed feedback ratchet 

The presence of a lag time in the control can cause negative currents and complicated dynamics that depends on 
the number of particles. Let us first study the few particle case (including one particle). 



1. Few particles 



When the control protocol presents a time delay the system performs worse because the delayed action of the 
controller implies some wrong actions. Moreover, for large time delays the controller is unable to surmount the 
potential shape asymmetry and eventually the net current becomes negative. See Figs. El and SI 

For increasing time delays the correlation between the present sign of the net force and the measured sign that the 
controller actually uses decreases. Thus the controller action begins to be uncorrelated to the present state of the 
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FIG. 2: Evolution of the net force per particle for the non delayed feedback ratchet for N — 20 (few particles) and A'^ = 10 
(many particles). The potential is "on" only in the time intervals such that /(f) is positive. Parameters of the potential: 
Vo = 5kBT and a = 2/3. Units: L = 1, D = 1, and ksT = 1. 




FIG. 3: Center-of-mass velocity (icm) as a function of the time delay r for different numbers A'^ of particles under the few 
particle regime. Parameters of the potential: Vo = lOfcsT and a — 2/3. Units: L — 1, D = 1, and ksT = 1. 

system and it effectively begins to act as an open-loop ratchet 112] ■ In fact, for large delays the correlation between 
the state of the system and the measured retarded state is negligible and the negative flux becomes independent of 
the delay, because the ratchet is effectively open-loop controlled; see Fig. [31 Therefore transport reversal appears here 
as a consequence of the competition between the asymmetry of the ratchet potential and the inherent asymmetry 
of the protocol. We stress that the influence of the asymmetry of the feedback protocol itself is not tuned here 
trivially, but changing the delay r in the control. Other ways of tuning the influence of the feedback protocol could 
not lead to current reversal. For example, it can be shown that a feedback protocol that switches on/off following 
the maximization protocol but with a probability of error < p < 1/2 does not enable negative fluxes even for 
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FIG. 4: One particle velocity versus the time delay for different heights Vo of the potential and asymmetry parameter a = 2/3. 
Units: L = 1, D = 1, and ksT = 1. 
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FIG. 5: Evolution of the net force per particle for time delays r — 0.01 (left) and r — 0.20 (right) in the many particle case 
(A'' = 10^). White background regions stand for "on" potential and gray background regions for "off" potential. Parameters 
of the potential: Vb = 5fcsT and a = 2/3. Units: i = 1, D = 1, and fc_gT — 1. 

asymmetries a > 1/2 (see Ref. Q). 

For the delayed protocol considered the critical value of the delay that gives zero current (thus the current is positive 
for smaller delays and negative for larger delays) is related with the characteristic time in which the information 
about the state of the system is effectively lost, so the delayed maximization protocol is not able to achieve its goal of 
producing a positive current for delays larger than the critical one. Increasing the height Vq of the potential implies 
a faster dynamics. Therefore the critical delay is expected to decrease with the height of the potential, in agreement 
with our simulations; see Fig. 3] It is important to note that this critical delay tends to a constant nonzero value as 
Vb — > oo. The reason is the same that makes the absolute value of the flux in both closed- loop and open-loop ratchets 
does not grow indefinitely as the potential goes up. These fluxes tend to a finite value because the time spent by the 
particles in diffusing during the off potential state goes to a constant in the absence of inertia. Note also that the 
critical delay decreases with the number of particles (see Fig. [31) 

Let us now study the many particle case, which exhibits a completely different dynamics. 



2. Many particles 



For large ensembles of particles {N > 10^ — 10'^) the flux is nearly zero in the non-delayed protocol (see Sec. 1111 A|) . 
but the introduction of a time delay stabilizes quasiperiodic solutions that give noticeable negative currents. We 
have found that, after a transient time, the delayed control allows the system to synchronize into a stable mode of 
oscillation such that the net force per particle evolves quasi-periodically. The evolution is not strictly periodic due to 
the stochastic nature of the dynamics. See Fig. [S] 

For small delays the on and off times of the non-delayed dynamics (Sec. 1111 Ap are enlarged owing to the delay, and 
the net force per particle evolves with a more regular pattern (Fig. [51 left). When the potential is switched on the 
net force per particle begins to diminish (because the potential asymmetry acts against the feedback protocol) and 
rapidly gets negative, but the potential still remains 'on' during a time r after the force changed its sign. On the 
other hand, when the potential is switched off the net force grows, becomes positive, and induces an 'on' switching a 
time r later. The result is a quasiperiodic dynamics with an small quasiperiod T > 2t; see Fig. [5l (left). We highhght 
that these types of solution are only observed for asymmetries a > 1/2, as they are consequences of the competing 
asymmetry of the potential; they do not appear for asymmetries a < 1/2 that support positive transport and exhibit 
a different behavior related with the enlargement of the tails of the net force per particle [12, [IjI ■ 

For larger delays there are stable solutions of quasiperiods T = 2r/(2ri-f 1), n = 0, 1, . . . , i.e., solutions that contain 
an odd number of semiperiods T/2 in the time delay r. See Fig. [5] (right) for instance. The competing asymmetry 
of the potential causes the stabilization of those solutions where, due to the delay, the controller switches on when 
the present net force is negative and switches off when it would be positive, that is, the controller acts contrary 
to its intentions and gives a negative flux. These branches are the counterparts of the solutions of quasiperiods 
T = r/n, n = 1,2, . . . (r containing an even number of semi-quasiperiods) observed for asymmetries a < 1/2 |12| . 
The difference of one semiperiod is due to the effectively reversed operation of the controller caused by the combined 
effect of the competing asymmetry of the potential and the delay. Some of these branches are plotted in Fig [HI for 
both the cases of cooperation (positive currents) and competition (negative currents) of asymmetries. 

It is important to note that the average velocity for all these branches can be reexpressed in terms of one of 
them. Let us define g{T) := (icm)r(''") as the average velocity for the branch of period T — t, which is present for 
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FIG. 6: Left panel: Center-of-mass velocity (icm) versus the delay t in the many particle case (A*' = IC). The region of small 
delays and the first two branches for asymmetries parameters a = 1/3 (positive flux) and a = 2/3 (negative flux) are plotted 
for height of the potential Vo — bknT. Units: L — 1 and ksT = 1. Right panel: First two branches for asymmetries a — 1/3 
and a = 2/3 for potential height Vo — bksT and A'' = 10^ particles, rescaled according to scaling laws © and (O- 

cooperative potential asymmetry. For these asymmetries, a < 1/2, we shovired in Ref. [12] that the average velocities 
of the branches of periods T = r/n are given by 



(ie,„>^(T)-.g(^) for a < 1/2. 



(5) 



On the contrary, for competing asymmetries, a > 1/2, we have found here that the solutions have quasiperiods 
T = 2T/(2ri + 1), and furthermore, the average velocities of these branches are given by 

(6) 



(icm) ,1^ (r) = -9 [^j for a > 1/2. 

Consequently, given one of the branches all the others can be predicted; see Fig. [5] This also implies that the analytical 
results obtained in Ref. [I4I for cooperative potential asymmetry (a < 1/2) are directly extended to the competing 
potential asymmetry case (a > 1/2), just using the relation in Eq. ([6]). 



IV. CONCLUSIONS 



We have studied the performance of feedback flashing ratchets when there is competition between the asymmetry in 
the potential and the asymmetry in the control protocol, and we have also studied the effects of tuning their relative 
influences in the dynamics by the introduction of a time delay. An experimental realization of a flashing ratchet 
has been performed in [16| | by using polystyrene latex spheres of diameters d ~ 0.25 — 1 fim in an aqueous solution 
[viscosity 77 ~ 10^'^ Pa ■ s; D = kBT/{2>Tir}dy\ exposed to a sawtooth dielectric potential of period i ~ 50 iim. This 
experimental setup can be modified to become an experimental realization of a feedback flashing ratchet by monitoring 
the particles with a conventional charge-coupled device (CCD) of about 30 fps and processing the images to switch on 
or off the ratchet potential in accordance with the particle positions. The time delays considered here are introduced 
by delaying the action of the controller a time between r = O.QIL^ / D ^ 10s and t — Q.5L^/D ~ 500s. (For a more 
detailed discussion see Ref. |l3l-) Indeed a sophisticated feedback control has been recently implemented in Ref. [1§ |. 
where images of a Brownian particle are acquired on a high-sensitivity CCD of up to 300 fps and thereafter a software 
processes the information to extract the position of the particle and apply a feedback voltage. On the other hand, 
we highlight that the viscous friction coefficient 7 depends on the shape and the size of the Brownian particle. Thus, 
as the adimensional delay must be multiplied by the factor L^/D = ^L'^/kgT in order to recover physical units, 
Brownian particles of different shape and size respond differently to a given time delay. This effect could be useful 
for separating different kinds of macromolecules. 

We have seen that the performance of the system with competing asymmetries differs significatively from its 
counterpart ratchet with cooperating asymmetries. In the absence of delay the competition of asymmetries implies a 
decay of the current with the size of the ensemble much stronger than in the cooperative case (l/\/]V vs 1/ IniV). In 
the presence of delay the dynamics becomes richer with a current reversal for large delays. In the few particle regime 
the change from positive to negative current can be understood as a change from a purely closed-loop control to an 
effective open-loop control. On the other hand, in the many particle case the negative current regime appears for 
large enough delays as the consequence of the stabilization of several branches of quasiperiodic solutions. These stable 
branches have the opposite sign and are one semiperiod displaced with respect to those obtained for cooperating 



asymmetries, they also have a direct relation with them that allows the extension for the competing asymmetries case 
of the analytical results found in Ref. ,12J for cooperative asymmetries. 
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